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Abstract 

We present the measurements of D°, D^, D*^ meson v'2 as well as D° 7?aa in different directions 
with respect to the estimated reaction plane in Pb-Pb collisions at yfs^ - 2.76 TeV with the 
ALICE detector at the LHC. 



1. Introduction 

On the basis of thermodynamic al considerations and Quantum Chromo-Dynamics (QCD) 
calculations, nuclear matter, in conditions of high temperature and density, is expected to undergo 
a phase transition to a deconfined state: the Quark Gluon Plasma (QGP) 1 1]. These conditions 
can be recreated in high energy nucleus-nucleus collisions. 

Non-central nucleus-nucleus collisions are characterized by an initial geometrical anisotropy 
with respect to the reaction plane (the plane defined by the beam direction and the impact pa- 
rameter). As a consequence of the different pressure gradients in the in-plane and out-of-plane 
regions, this spatial anisotropy is converted to a momentum anisotropy of produced 1ow-/?t par- 
ticles. The momentum anisotropy can be quantified via a Fourier expansion of the azimuthal 
angle with respect to the estimated reaction plane. The second coefficient of this expansion V2 is 
called elliptic flow [2] . The azimuthal anisotropy of heavy-flavour hadron production is sensitive 
to the degree of thermalization of heavy quarks in the expanding medium and to the path length 
dependence of their energy loss iSH, S la 0, Si ■ 

2. Data sample and event plane determination 

The D meson azimuthal anisotropy measurements presented in this report are performed with 
the ALICE detector using the 201 1 data sample collected with a minimum bias trigger given 
by the coincidence of signals in the VZERO scintillators and the Silicon Pixel Detector (SPD) 
and a trigger tuned to enhance the sample of 50% most central collisions, based on the VZERO 
signal amplitude, which is used to classify the events according to centrality. The analyses have 
been peroformed with 9.5 x 10^ Pb-Pb collisions in the centrality class 30-50%, 7.1 x 10^ in 
15-30% and 16x10'' in 0-7.5%. 

The azimuthal direction of the event plane, which estimates the reaction plane direction, is 
determinated from the distribution of charged tracks in the Q < rf < 0.8 pseudo-rapidity interval 
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of the Time Projection Chamber (TPC) using the equation 
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where *P is the second harmonic event plane and (pi is the angle of the track in the ALICE 
reference frame. The event plane resolution {R2) was computed considering each event as splitted 
in 2 sub-events, made of randomly combined tracks and associated to the correspondent event 
plane angle, following the prescription in liotl . The resulting event plane is flat to about 2% level 
and its resolution is 0.86 in the 30-50% centrality class, 0.9 in 15-30% and 0.75 in 0-7.5%. 

3. D meson reconstruction, V2 and /?aa determination 

D meson reconstruction in ALICE is based on the invariant mass analysis of fully recon- 
structed hadronic decay topologies of these channels: D" — » K ;7r+, — > K ;r^7r^ and D*^ — » D'^n* 
The separation of about a few hundred /xm between the primary and the secondary vertices, pe- 
culiar of D*^ and decay, is exploited to reduce the combinatorial background. In ALICE, 
charged tracks are reconstructed with a six layers silicon tracker (Inner Tracking System) and 
the Time Projection Chamber (TPC). In particular, the Silicon Pixel Detector (SPD) provides a 
measurement of track impact parameter to the primary vertex with a resolution better than 65 jjm 
for tracks with pj > 1 GeV/c. Particle identification is provided by the time of flight measure- 
ment in the Time-Of-Flight detector for kaons with pr < 2 GeV/c and by the specific energy 
deposit in the TPC. For high-px particles, a selection is applied in order to reject protons with- 
out loosing possible kaons. The contribution of signal candidates with wrong mass assignement 
to the final state hadrons was found to be negligible and it does not bias the signal extraction. 
Further details on the analysis can be found in ill ill . The invariant mass distribution of recon- 
structed D meson candidates is considered in the in-plane and out-of-plane regions defined as 
[0 < < 1) U < A(/? < tt) and [| < A(/j < ^) respectively, where Aip is the azimuthal angle 
of the reconstructed candidate with respect to the event plane 'F. An invariant mass analysis, 
based on a Gaussian fit, is used to obtain the signal yield in the two regions (A^in, A^out)- Assum- 
ing that the reconstruction and selection efficiency is indipendent of Ay?, it is possible to directly 
compute the elliptic flow from the formula (jH left), where R2 is the event plane resolution. 

n 1 Mn-A^out „ , , dNAA/dpj 
^2 - TTTTT Tt ' '^l^AKPl) - — — : — — (^) 

AR2Nm+NouY (Taa) x do-pp/dp-r 

The /?AA measurement is based on the comparison between the yield measured in AA collisions 
and the cross section measured in pp collisions, scaled by the overlap nuclear function, as shown 
in formula (|2] right). For this analysis, the yields A^in and A^out were corrected for acceptance 
and efficiency as a function of pj, using PYTHIA -1- HIJING Monte Carlo simulations il2ll . 
Details on the correction for the secondary D meson and the pp measurement, used as reference, 
can be found in [12]. 

Several sources of systematic uncertainties were considered in both analyses. The main con- 
tributions come from uncertainties on the yield extraction and from the topological cut selection. 
The former is estimated using different background functions in the fitting precedure, different 
mass range, and a bin counting method. The analysis was also repeated with three different sets 
of topological cuts. Yield extraction and cut variation systematic uncertainties are estimated to 
be at maximum 0.05 absolute value, for the V2 analysis. The signal sample considered contains 
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Figure 1: (left) D", D+ and D*^ i'2 measured by ALICE in the 30-50% centrality class. Empty boxes show systematic 
uncertaintites from data and the shaded area the B feed-one one. Black points shows the charged hadron V2. (light) 
D" V2 as a function of centrality for the bins 2 < pT < 3 GeV/c, 3 < pj < 4 GeV/c. Empty boxes show systematic 
uncertaintites from data and the shaded areas the B feed-down one. 

a fraction of D mesons coming from B decays, thus the measured elliptic flow is a combination 
of prompt and secondary D meson anisotropy. The feed-down correction uses as input pQCD 
calculations of B production, as explained in L 1 lil . The resulting systematic uncertainty on V2 is 
up to +23%, including a conservative variation of the unknown R\a and V2 of D meson from B 
decays. The centrality dependence of the event plane resolution in the centrality range consid- 
ered was also taken into account, together with the uncertainty on the resolution. For the 30-50% 
centrality class, these uncertainties were found to be ±3% and +7% respectively. The small dif- 
ference of efficiencies in the two A(f region was found to give a negligible contribution to the V2 
uncertainties. 



The D meson anisotropy was measured in the 30-50% centrality class: for D in the range 2 < 
Pt < 16 GeV/c, D+ in 3 < /?T < 8 GeV/c and D*+ in 2 < /?t < 20 GeV/c. The three results are 
compatible within statistical uncertainties and all are compatible with the ALICE measurement 
of charged hadrons in the same rapidity region (Fig. [1] left). The D" V2 measurement was also 
performed in the centrality classes 15-30% and 0-7.5% in the range 2 < pj < 16 GeV/c. The 
results for the first two pj intervals show a hint of increasing V2 from central to semi-peripheral 
collisions, as reported in Fig[T] right. V2 was measured using two-particles correlation methods 
as well. The results are consistent with those based on the event plane. The measurement of 
the D meson anisotropy was compared to theoretical predictions il3ll : it was observed that a 
simultaneous description of the anisotropy and of the Raa suppression is challenging for models. 

The D" Raa in the two azimuthal regions, in-plane and out-of-plane, was measured in the 
range 2 < pj < 16 GeV/c. The measurement shows more suppression in the out-of-plane region 
with respect to the in-plane one, suggesting a path length dependence of heavy quark energy loss, 
which is expected to be the dominant effect at high pj (Fig.|2]). At low pj, elliptic flow may also 
contribute to the observed anisotropy. In Fig. [2] we compare the D" /?aa in the two azimuthal 
regions with theoretical predictions. The WHDG |14] and POWLANG |5] models show a good 
agreement with out-of-plane results while they somewhat oversuppress the in-plane yields. This 



4. Results 
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Figure 2: D" Raa vs event plane for the 30-50% centrality class. Empty boxes show the uncoiTelated systematic uncer- 
tainties between the two measurements, the empty brackets the con'elated systematic uncertainties that would shift both 
measurements in the same direction and shaded areas the anticorrelated uncertainties that would shift the measurements 
in opposite directions. 



high-/7T feature is present also for the BAMPS model which instead can describe the data in 
the low-pj range. 

5. Conclusions 

The measurement of D meson V2 in Pb-Pb collisions at -y/iNN = 2.76 TeV with ALICE 
was presented. The results indicate V2 > in the pj range 2 < pj < 6 GeV/c with a 3cr 
significance in the 30-50% centrality class. A hint of centrality dependence is observed. The 
measurement of the D" R^a around two orthogonal directions with respect to the event plane, in 
the 30-50% centrality class, indicates a larger suppression in the out-of-plane direction, where 
the initial overlap region of the two nuclei is more elonged, suggesting a path lenght dependence 
of heavy-quark energy loss. 
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